Fermi hypernetted-chain study of half-filled Landau levels with broken rotational 

symmetry 



o 
o 

(N 

> 
O 



i 

c 

o 
o 



> 
o 

(N 



o3 



i 

c 

o 
o 



Orion Ciftja and Carlos Wexler 
Department of Physics and Astronomy, University of Missouri-Columbia, Columbia, Missouri 65211, USA 

(February 1, 2008) 

We investigate broken rotational symmetry (BRS) states at half-filling of the valence Landau 
level (LL). We generalize Rezayi- Read's (RR) composite fermion (CF) trial wavefunctions to include 
anisotropic coupling of the flux quanta to electrons, thus generating a nematic order in the underlying 
CF liquid. Using the Fermi hypernetted-chain (FHNC) method which readily gives results in the 
thermodynamic limit, we determine in detail the properties of these states. By using the anisotropic 
pair distribution and static structure functions we determine the correlation energy and find that, 
as expected, RR's state is stable in the lowest LL, whereas BRS states may occur at half-filling of 
higher LLs, with possible connection to the recently discovered quantum Hall liquid crystals. 

PACS: 73.43.-f, 73.20.Dx, 64.70.Md. 



I. INTRODUCTION 

During the past two decades the physics of two- 
dimensional electron systems (2DES) has provided some 
of the richest grounds for new developments in condensed 
matter physics. In particular the integer El and frac- 
tional Q quantum Hall effects (QH) in strong magnetic 
fields are some of the most remarkable phenomena dis- 
covered in the second half of the XX Century and rival 
superconductivity in their fundamental significance by 
manifesting quantum mechanics on macroscopic scales 
and providing a major impetus to the development of 
new ideas in many-body physics H , like the existence of 
fractionally charged quasiparticles topological quan- 
tum numbers Q, chiral Luttinger liquids q|, composite 
particles j7j , etc. One of the reasons that 2DES keep sup- 
plying novel and exciting results is the improved quality 
of the samples with mobilities increasing roughly expo- 
nentially with time, thus allowing the emergence of new, 
subtler effects due to electronic correlations (which are 
enhanced because of the reduced dimensionality). 

Consider initially a 2DES at half-filling (y = 1/2) of 
the lowest Landau level (LLL) . Experimentally, this state 
does not exhibit the typical features of the fractional 
quantum Hall effect (FQHE), namely the very precise 
quantization of the transverse conductance in units of 
e 2 /h, or the vanishing longitudinal resistivity. Yet, the 
resistivity shows a broad minimum ||, and anomalous 
behavior in the propagation of surface acoustic waves 
H has been observed. Early numerical work by Hal- 
dane Jl(j], suggested that v = 1/2 is not incompress- 
ible. The overall experimental evidence strongly suggests 
that, in the LLL, for half- filling the system behaves like 
a strange Fermi liquid at close to zero effective magnetic 
field 

A theory of compressible Fermi-liquid- like behavior at 
half filling has been proposed by Halperin, Lee, and Read 
(HLR) |Q : a two-dimensional (2D) system of electrons 
subjected to an external perpendicular magnetic field, at 
half filling of the LLL can be transformed to a mathemat- 
ically equivalent system of fcrmions interacting with a 



Chern-Simons gauge field such that the average effective 
magnetic field acting on the fermions is zero [ pd[ . Since 
these fermions do not see a net magnetic field, they can 
form a 2D Fermi sea of uniform density. In a very suc- 
cessful approach, based on the composite fermion (CF) 
theory fi|i5|], Rezayi- Read |16| (RR) described the half- 
filled state by a correlated Fermi wave function that is 
a product of a Slater determinant of plane waves with a 
Jastrow factor corresponding to a Bose Laughlin state at 
half filling: 
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where </5k(i"i) are 2D plane waves for fully spin polarized 
CF states that fill a 2D disk in reciprocal space with 
Fermi momentum hp, Zj = Xj + iyj is the complex 2D 
coordinate of j-th electron and Po is a projector onto the 
LLL (L = 0). We work in units of the magnetic length 
{ll = H/eB = 1). 

The situation is dramatically different for the next LL 
(LL index L = 1), where for v ~ 5/2 a quantized Hall 
conductance is observed along with a strong reduction of 
the longitudinal conductance [p"7|] , which are an indica- 
tion of a Cooper pairing instability of the CFs A 
discussion of this state is a wide and complex topic by 
itself and this, along with the properties of the v = 3/2 
and 7/2 states (half filling of the upper spin sub-band) 
will not be addressed in this work. 

In this paper we discuss some aspects related to the 
many interesting new phenomena that have recently 
emerged in nearly half-filled higher LLs (with LL index 
L > 2), in particular the extreme anisotropy measured 
in the low temperature magnetotransport [p_Q| — 121| . This 
anisotropic behavior has been attributed to the formation 
of a nematic phase of the 2DES which, at higher tempera- 
tures undergoes a nematic to isotropic transition p2],p3[. 
The motivation of our work is to study these nematic 
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phases by means of many-body trial wavefunctions with 
broken rotation symmetry (BRS) at half-filling of a Lan- 
dau level. Previously H] we reported on the existence of 
a BRS instability of the Laughlin state M for 1/3-fillcd 
higher LLs (with L > 1). Here we discuss a similar pro- 
cedure for the more complex case of 1/2-filling, which is 
more closely related to the experiments on anisotropic 



phases fl9| |2l| 



Similarly to the Musaelian and Joynt's (MJ) [|25[ gen- 
eralization of Laughlin's FQHE state Q used by us in 
Rcf. 1 24 1, we add a symmetry breaking parameter a in 
the RR wavefunction [Eq. (pi)] : 



A' 
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where Pj, now represents a projector onto the L 
LL. This wavefunction represents a homogeneous Fermi- 
liquid-like state at half filling factor, and for a / has 
nematic order (for a = we recover the RR wavefunction 
which is obviously isotropic) . Note that the magnitude of 
a is related to the amount of anisotropy and its phase to 
the angle the director of the nematic (for real a the sys- 
tem will have a stronger modulation in the x-direction, 
and therefore likely have larger conductance in the per- 
pendicular direction: a yy > o~ xx ). This wavefunction is 
an obvious starting point to study the nematic quantum 
Hall liquid crystals at half filling, by facilitating the sys- 
tematic study of the energy dependence of BRS states 
for diverse physical parameters (LL index, width of the 
2DES, etc). 

We study the BRS state for 1/2-filling of the valence 
LL (i.e. v = M + 1/2 with M integer) by using the 
Fermi hypernetted-chain (FHNC) method This 
method allows us to compute physical quantities in the 
thermodynamic limit, without the limitations of using a 
finite number of particles that hinder other techniques, 
where the extrapolation of the results to the thermody- 
namic limit is not totally unambiguous. 

We find that, for realistic interaction potentials, the 
RR (isotropic) state is stable in the LLL, whereas a BRS 
state is possible in higher LLs. While this is consistent 
with the view of BRS states for L > 2, it does not re- 
flect the situation for L = 1, where it is believed that 
CF-s can form an incompressible state by a Bardeen- 
Cooper-Schrieffer (BCS)-like pairing J3l| as first pointed 
out by Moore and Read |3^j . Recent exact diagonaliza- 
tion results by Morf J|3| strongly suggest that such spin- 
polarized Pfaffian state is the best candidate to describe 
this filling. 

In Sec. H we present the basic theoretical calculations 
needed to determine the stability of an isotropic or BRS 
state. A detailed description of the FHNC formalism in 
the context of the BRS wave function [Eq. (||)] is given 
in Sec. III. The results for the BRS state in the LLL and 



their extension in higher LLs are discussed in Sec. IV 



II. BASIC THEORY 

In this work we study the stability of different states 
by using trial wavefunctions of the form of Eq. @. We 
perform this analysis by comparing the energy in each of 
these states to find the optimum value for the anisotropy- 
generating parameter a. The potential, or correlation 
energy per electron is given by: 



1 (go|Vlgq) 



£ / d 2 rV(r) [g(r)-l] , (3) 



where V represents the electron-electron, electron- 
background, and background-background interaction; 
and g(r) is the (angle-dependent) pair distribution func- 
tion given by 



ff( r ) 
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where r = V2 — v\ . The following normalization con- 
dition, p J d 2 r [g(r) — 1] = —1, can be used as a con- 
venient check for numerical procedures. For an ideal 
2D sample the interaction is a pure Coulomb potential 
V(r) ~ e 2 /(er) , while in samples with finite thickness a 
reasonable choice is the Zhang Das Sarma (ZDS) poten- 
tial pi V(r) = e 2 /(e v / r 2 + A 2 ), where A is of the order 
of the sample thickness. Alternatively, the correlation 
energy can be computed in reciprocal space: 



d 2 q ~ 
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V(q) [5(q) - 1] 



(5) 



where V(q) is the 2D Fourier transform Q (FT) of V(r) 
and 5(q) is the static structure factor: 



S(cO - 1 = P FT[g(r) - 1} 



(6) 



While both g(r) and 5(q) are angle-dependent (e.g. 
see Figs, [l] and ||), because the interaction potential is 
centrally symmetric, the energy E a depends only on the 
angle-averaged pair distribution function or static struc- 
ture factor defined as: 
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The determination of either the pair distribution func- 
tion or the structure factor is generally a complicated 
integral problem that needs to be solved for each LL. 
However, its is known that if transitions to other LLs 
are neglected (i.e. a single-LL approximation), g(r) and 
5(q) at higher LL are simply related to those at the LLL 
(L = 0) by means of a convolution or product respec- 
tively. We will apply this approximation (which, more- 
over, quenches the kinetic energy in higher LLs as well). 
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It is then, sufficient to compute these distribution func- 
tions once in the LLL and then the correlation energy 
per electron is given by 



If, 



d 2 q - 
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V cS (q) [S(q) - 1] 
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where V B g(q) = V(q) [L ^(q 2 / '2)\ 2 . Ll(z) are Laguerre 
polynomials, and S(q) is calculated in the LLL (L = 0). 

In what follows we compute <7(r) and S(q) using the 
FHNC method. 



III. THE FERMI HNC METHOD FOR THE 
BROKEN ROTATIONAL SYMMETRY STATE 

The development of the FHNC method for Fermi sys- 
tems |26|,^7j allows one to estimate accurately the expec- 
tation value of a Hamiltonian, the pair distribution func- 
tion and related quantities associated with a Jastrow- 
Slater wave function and other more complex many-body 
Fermi wave functions. The FHNC method treats the 
correlated system of particles a priori in the thermody- 
namic limit therefore is extremely useful on the study 
of infinite homogeneous Fermi systems described by cor- 
related Fermi wave functions, having (but not limited 
to) a Jastrow-Slater form. In addition one can prove 
that the FHNC scheme achieves convergence of the ex- 
pectation value of a Hamiltonian within some expansion 
scheme j3f| thus substantiating the belief that FHNC 
gives an accurate upper bound for the energy and other 
related quantities pq,[3j|. 

An important problem arises at this stage. The pro- 
jection onto the L th LL performed in the RR wavefunc- 



y] leads to a wavefunction 



tion or its generalization [Eq. (^ 
that cannot be directly treated within the FHNC formal- 
ism, because the simple Jastrow x Slater determinant 
structure of single particle orbitals is lost. We therefore 
use an unprojected version of Eq. (|J), which although 
approximate, is believed to contain the most important 
physics, especially since the Jastrow factors already sig- 
nificantly annihilate higher-LL components of the wave- 
function pq| . In addition, although the wave function 
given in Eq. (^) has a Jastrow-Slater form after dropping 
the projector Pl, the FHNC method in this case differs 
from the standard approach since two-body correlation 
factors and related quantities depend not only on inter- 
particle distance, but also on the relative angle between 
the particles. In order to calculate the pair distribution 
function the modulus square of the wave function: 

|* Q (r 1 ,...,r iv )| 2 = e&«^-^ ( 9 ) 
x e -£, =1 ^Hdet fckfa)] | 2 , 

where u(z) = In \z — a\ 2 + In \z + a\ 2 , is expanded pertur- 
batively in terms of h(rij) = exp[u(r,y)] — 1 and may be 
ordered as: 



N N N 
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x|det[^ k (r,)] | 2 . 



(10) 



In addition to dynamical correlations between particles, 
described by the pseudopotential u(zi-Zj), there are also 
statistical correlations described by the Slater determi- 
nant which renders the whole state antisymmetric. Sim- 
ilar to the Jastrow part, we may expand |det [^k(**)] | 2 
in the number of exchanges between particles. In this 
way Eq. ( [l0|) becomes a quite symmetrical expansion in 
the number of dynamical correlation factors, /i(r^), as 
well as the number of statistical correlation factors. The 
resulting cluster terms in the expression for g(r±2) con- 
tain both kinds of correlations and may be represented by 
cluster diagrams. As in the Bose case, the associated pair 
distribution function <?(i"i2) is then given by the sum of 
all linked irreducible diagrams obeying well-defined topo- 
logical rules p6| . 

One defines nodal, composite (non-nodal), and elemen- 
tary diagrams as in the Bose case, but there are now 
four different types for each of them. The four differ- 
ent classes of nodal, composite and elementary diagrams 
are generally denoted by dd (direct-direct), de (direct- 
exchange), ee ( exchange- exchange ) , and cc (circular- 
exchange). Then the pair distribution function is ob- 
tained from the set of FHNC equations given in the Ap- 
pendix, where p is the density, g s {— 1) is the spin degen- 
eracy and l{r 12) is the familiar statistical exchange factor 
for the 2D Slater determinant given by 



/(ri,r 2 ) =2 



■hjkF r 12 ) 
kF ri2 



(11) 



where r 12 = |r 2 — r x | and J\(x) is the first order Bessel 
function. 

For standard systems like the 3 He Fermi liquid, the 
pair correlation factor is short-range and heals to 1 for 
large distances, therefore the function exp[u(r^)] — 1 
provides a possible expansion parameter [note that in 
order to apply the Fermi HNC expansion, the corre- 
lation (pseudo) potential has to satisfy the conditions: 
u(ry — * 0) — » —00 and u(ry — * +00) — > 0]. In the 
case of the BRS wave function, the (pseudo) potential 
u(v) is logarithmically long-range, however formally it is 
possible to extend the method by splitting the pseudopo- 
tential associated with the Jastrow part into a short- and 
a long-ranged part, respectively, 

u(r 12 ) = u s (r 12 ) + ui(r 12 ) , (12) 

with the dd nodal and composite function similarly split: 

N M {r 12 ) = N dds (r 12 ) - ui(r l2 ) , (13) 
X dd (r 12 ) = X dds {r 12 ) + u t {r 12 ) . (14) 

The splitting is done subject to the conditions: 
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t»(ria) + N dd (r 12 ) = u s (r 12 ) + N dds (r 12 ) , (15) 
N dd (r 12 ) + X dd (r 12 ) = N dds (r l2 ) + X dds (r 12 ) . (16) 

The short-range function U s (ri 2 ) (going to -oo for small 
distances and healing to for large distances) and its its 
long-range counterpart are then chosen as: 

u,(r 12 ) = -2K (Q\r 12 -a\) - 2 K (Q\r 12 + a\) , (17) 
ui(r 12 ) = 2 [ln(|r 12 - a\) + K (Q\r 12 - a\)} 

+2 [ln(|r 12 + a\) + K {Q\r 12 + a\)} , (18) 

where Kq(x) is the modified Bcssel function, and Q is 
a cut-off parameter of order 1. We recall that the 2D 
FT f§ of ui{r 12 ) is: 



ui(q) = - 2 / 2 
q 2 {q 2 



4ttQ 2 



Q 2 ) 



(19) 



The final set of equations is solved by a standard iter- 
ative procedure. There is one necessary approximation 
within the FHNC method in order to obtain a closed 
set of equations for the nodal and non-nodal functions: 
a small set of cluster diagrams (corresponding to the 
so called elementary diagrams), which cannot be in- 
cluded ab initio in the method, needs to be somehow 
estimated outside the FHNC method. Several schemes 
have been devised to include the contribution of such di- 
agrams at various levels of approximation, however the 
simplest approximation of totally neglecting these terms 
(called the FHNC/0 approximation where we assume 
E dd = E de = E ee = E cc = 0) generally leads to very 
reliable results |2q, ^6[^,|3^] and we have adopted it in 
this paper. 

IV. RESULTS AND DISCUSSIONS 

In this work we applied the FHNC theory to study 
the BRS state at filling 1/2 of an arbitrary LL. For the 
sake of simplicity we neglected the elementary diagrams 
(FHNC/0). This has allowed us to determine to a rea- 
sonable accuracy the pair distribution function and the 
static structure factor. In order to compare the a = 
(RR state, isotropic) with the a ^ (BRS state, or ne- 
matic) we studied the properties of the BRS wave func- 
tion for several as with magnitudes between and 3 (in 
general a = |a|e iea , but without loosing generality we 
considered only 8 a = 0). 

In Fig. ^we plot the pair distribution function g (r) for 
a = 2 [panels (a) and (b)], and the angle-averaged pair 
distribution function g(r) corresponding to a = 0, 1, 2, 
2.5 and 3 [panels (c) and (d)]. It is interesting to note, 
for a / 0, the noticeable angle-dependence of g(r), and 
the splitting of the triple node at the origin to a simple 
node at the origin and additional simple nodes at r = a 
and angle 8 = 6 a , 8 a + ir (6 a = in this case). From 
the energetics point of view [see Eq. (g)], a major con- 
sideration is the strong dependence of g(r) on the value 
of the parameter a. As a is increased the major peak 



of g~(r) becomes less pronounced and for a w 3 it devel- 
ops a shoulder for small r [panel (c)]. In addition, note 
change in the small-r behavior of g(r) which switches 



(for a = 0) to 



as a is increased [panel 



from o 

(d)]. In general, for small-r, g(r) has almost no angular 
dependence jl0|, and for a ^ 0, g(r w 0, 6) ~ C a r 2 for 
< r < 0.5, where C a ~ 0.024 a 17 . 

In Fig. H we plot the static structure factor 5(q) for 
a — 2 (top panel), where the most important feature is 
the emergence of peaks in S(q) characteristic of a ne- 
matic structure; and the angle-averaged static structure 
factor S(q) corresponding to a — 0, 1, 2, 2.5 and 3 (bot- 
tom panel). Note the considerable dependence of S(q) on 
a: as it increases the peak is broadened and flattened, 
with no significant change in the small-g behavior. 

One can compute the correlation energy per particle 
either directly from Eqs. (§),(§), or (||) to determine the 
energy per electron for arbitrary values of the Fermi BRS 
parameter a, the 2D system width A, and Landau level 
index L. The following simplified formula can be used in 
view of Eq. (0) : 



E L a {\) 



1 

4tt 



dqqV(q,X) [L L {\)] 2 [S{q) 1] 



(20) 



where V(q, A) = (27re 2 /eg) exp(-Ag) is the 2D FT of 
the ZDS interaction potential |34j. In addition to allow- 
ing straightforward calculations to be extended to any 
LL, Eq. ( pp| ) permits higher numerical accuracy on the 
calculation of E a since S(q) saturates exponentially to 1 
for relatively small values of q as compared to ~g(r) . 

Fig. H shows the energy difference between BRS states 
with a — 1, 2, 2.5 and 3, and the isotropic state with 
a = 0. Our findings indicate that in the LLL (L = 0) 
the Rezayi-Read (RR) state is stable for any A, since all 
a^O states have higher energies (top panel). 

The situation changes considerably in higher LLs (L > 
1). In general, Fermi BRS states are found to have lower 
energies and the compressible RR state is unstable to- 
wards a nematic state (see lower panels of Fig. ||) . Con- 
trary to our findings for 1 /3-fillcd LLs |24| , we do not see 
a runaway instability but rather find that there are opti- 
mal values: a\ w 2 for the first LL (L = 1), and a 2 « 1 
for L = 2. Whereas no BRS state has been so far ob- 
served for L = 1 (likely because another correlated state 
not considered here is lower in energy) it is likely that the 
BRS state found here is related to the low-temperature 
anisotropic conductance found in higher LLs [|19| - |23| . 

At this point it is important to comment on how precise 
our determination of these energy differences is. There 
are two aspects to this problem. First, the reader should 
note that we used an unprotected wave function in place 
of the RR generalization [Eq. (||)]. While it would be 
highly desirable to incorporate this LL projection oper- 
ators into the formalism this is not feasible within the 
FHNC method used here. We should note, however, that 
the presence of the Jastrow factors in Eq. (||), already 
provides a considerable projection into the LLL, and is 
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believed to be particularly effective as far as ground-state 
properties are concerned |Q. The second aspect is that 
the FHNC /0 method is essentially a variational method, 
giving energies that constitute an upper bound to the ex- 
act ground state energy j4l[] . While a precise estimate of 
these errors can only be made a posteriori by comparing 
these results with alternative calculations, these methods 
have proven to be quite reliable in scenarios ranging from 
other partially filled LLs to 3 He systems 

While the energy differences we are interested (see Fig. 
^) are quite small, and perhaps smaller than the errors 
in the absolute values of the correlation energy in the 
states, we remark that these are not uncorrelated errors 
but systematic deviations due to the nature of the ap- 
proximations used, and energy differences will likely be 
considerably more precise. 

In conclusion, we applied the FHNC theory to study 
possible Fermi BRS states in a 1/2-filled LL. We find 
that the isotropic RR state is stable in the LLL for re- 
alistic interaction potentials. In higher LLs, BRS states 
with nematic order are energetically more favorable than 
the RR state, perhaps with a direct connection to the 
anisotropic states observed recently in high LLs 

We would like to acknowledge helpful discussions with 
A.T. Dorsey, M. Fogler, L. Radzihovsky and G. Vignale. 
This work was supported by the University of Missouri 
Research Board. 

APPENDIX A: 

For a Fermi system at density p and spin degeneracy 
g s (1 or 2) the sum of non- nodal (composite) diagrams is 
given by: 

X dd (r 12 ) = e ^)+N dd (r 12 )+E dd (r 12 ) _ Ndd (v 12 ) - 1 

(Al) 

X de (v 12 ) — e u ( ri2 )+JVdd(ri2)+-B £i d(ri2) 

x [A^ e (r 12 ) + E de (r 12 )} - N de {r 12 ) , (A2) 

X ee (v 12 ) = e «(r 12 )+^(r 13 )+^ d (r 12 )[ iVee ( ri2 ) 

+E ee (r 12 ) + \N de (r 12 ) + E de (r 12 )\ 2 - g s \N cc (r 12 ) 
+E cc {v 12 ) - l(r 12 )/g s \ 2 } ~ N ee (r 12 ) , (A3) 

A cc (r 12 ) = e ^i2) + N dd (r 12 ) + E dd (r 12 ) [Ncc{ri2) + 

-l{r 12 )/g s ] + l{v 12 )/g s - N cc {r 12 ) . (A4) 

The chain formation of the nodal diagrams is generated 
by convolution equations 

N dd (r 12 ) = p J dr 3 [X dd (r 13 ) + N dd {r 13 )}P(r 32 ) , (A5) 

N de (r 12 ) =p J dr 3 [X dd (r 13 )X ee {r 32 ) (A6) 
-X de (r 13 )X de (r 32 ) + [X de (r 13 ) + JV«te(r 13 )]P(r32)], 



N ee (r 12 ) = P J dr 3 [X de (r 13 )X de (r 32 ) (A7) 
-X dd {r l3 )X ee {r 32 ) + [X ee (r 13 ) + N ee {r 13 )]P{r 32 )} , 

N cc {r 12 ) = p J dr 3 [~l(r 13 )/g s + X cc (r 13 ) 

+N cc (r 13 )]X cc (r 32 ) , (A8) 

where 

P(r« ) - X dd {Tij ) + 2X de (r« ) ( A9) 

+p / dr k [X dd (Y ik )X ee (v kj ) - X de (r ik )X de (r kj )]. 



The pair distribution function is then given from: 

g(r 12 ) = 1 + X dd (r 12 ) + N dd (rl2) (A10) 
+2 [X de (r 12 ) + N de (r 12 )} + X ee (r 12 ) + N ee {r 12 ) . 
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FIG. 1. Pair distribution function for the BRS state at v = 1/2. (a) a = 2, surface plot of g(r) (the surface for y < was 
removed for clarity); (b) a — 2, dotted lines: g(r,8) for various 9 £ [0, 27r], full line: angle averaged g(r); (cj Angle averaged 
g(r) for a = 0, 1, 2, 2.5 and 3; fdj Small-r behavior of g(r), lines are fitting curves. Note the discrete nodes of g(r,0) at r = a, 
9 — 9 a ,6 a + -K (9 a — in this case). Calculations were performed in the FHNC/0 approximation. 
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FIG. 2. Static structure factor for the BRS state at 
v — 1/2. (a) a = 2, surface plot of 5*(q) (the surface for 
q v < was removed for clarity); (b) Angle averaged S(q) for 
a = 0, 1, 2, 2.5 and 3. Note the presence of peaks in S'(q) 
consistent with a nematic structure. Calculations were per- 
formed in the FHNC/0 approximation. 



FIG. 3. Energy per particle of the Fermi BRS states with 
a = 1, 2, 2.5 and 3 relative to the isotropic (a — 0) state: 
AE a (X) = E a (\)-E (\) for the LLL (L = 0) and higher LL-s 
as function of the short distance cut-off parameter A. Ener- 
gies are in units of e 2 /(d ). Note that in the LLL, the Fermi 
BRS states always have an energy higher than the isotropic 
state, whereas in higher LLs (L = 1, 2) there are ranges of A 
for which Fermi BRS states are favorable. 
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